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Abstract: Encarsia Forster consists of important parasitoids of whitefly (Bemisia tabaci) pests, including 
E. bimaculata, E. formosa and E. sophia, the three most important aphelinid parasitoids in China. 
Eight populations of Encarsia from the South, Southeast, North and Southwest of China, as well as two 
populations from Malaysia and Egypt, respectively , were collected in the present study, and their inter- 
species phylogenetic relationships were analyzed based on 28S rRNA D2 and D3 expansion regions. The 
D2 and D3 regions were consistent with each other, confirmed a closer genetic relationship between E. 
sophia and E. bimaculata since they both belong to the Encarisa strenus species group, compared to those 
between these two species and En. formosa. Results of the genetic distance analysis using 28S rRNA D2 
sequences revealed that there are certain genetic divergences within single species of the Encarsia 
parasitoids. The Guangzhou population of Encarsia sophia is more close to populations from Australia, 
Spain, Egypt and Ethiopia, but further from the population from Thailand. E. bimaculata populations 
from Sudan, Egypt and Guatemala as well as one population from Australia cluster together, while E. 
formosa Hengshui and Kunming populations cluster together with those from USA, UK and Greece, but 
are further from the Egypt population. The reasons for the inconsistency between the genetic and 


geographical distances of the Encarsia species are discussed. 
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1 INTRODUCTION 


The sweetpotato whitefly Bemisia tabaci 
(Hemiptera; Aleyrodidae ) is a globally destructive 
agricultural pest (De Barro, 1995; Liu et al., 2007; 
De Barro et al., 2011). The first record of B. tabaci 
in China was in the late 1940’ s (Chou, 1949) , but 
it was not considered as a serious pest until outbreaks 
in the mid-1990’ s due to the invasion of biotype B 
belonging to the Middle East-Asia Minor 1 genetic 
group (Qiu et al., 2003, 2007a; Wu et al., 2003 ) 
and then the subsequent invasion of the Q biotype 
belonging to the Mediterranean genetic group ( Chu 
et al., 2006). Since then, B. tabaci has caused 


serious damage to a wide range of vegetables, fiber 
and ornamental crops in more than 20 provinces 


across China ( Zhang et al., 2005; Li et al., 2011). 





One of the main challenges in B. tabaci 
management is that this pest continues to show an 
increasing level of resistance to a range of chemical 
acetamiprid , 


insecticides, including nitenpyram, 


imidacloprid and several pyrethroids ( Morin et al., 
2002; Horowitz et al., 2005; Byrne et al., 2010; 
Feng et al., 2010; Schuster et al., 2010). Several 
recent studies in China have shown that most of the 
routine pesticides available have lost their 
effectiveness in B. tabaci control (He et al., 2009; 
Qiu et al., 2009; Wang et al., 2009). Therefore, 
long-term sustainable management of B. tabaci 
requires an integrated approach in which biological 
control with natural enemies is one of the key 
components ( De Barro et al., 2000; Cuthbertson et 
al., 2011). 

During the last two decades, biological control 
agents of B. tabaci, including parasitoids, predators 
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and both entomopathogenic fungi and nematodes as 
well as their associated technologies, have become 
the focus of research in B. tabaci management 
studies ( Gerling et al., 2001; Ren et al., 2001; 
Cuthbertson et al., 2008, 2011; Stansly and 
Naranjo, 2010). In China, 54 species of aphelinid 
parasitoids of B. tabaci have been recorded in a 
series of extensive field surveys, including 41 
species of Encarsia and 12 species of Eretmocerus 
(Li et al., 2011). Three species of Encarsia, for 
E. formosa and E. 


sophia, are the dominant species in China as 


example, E. bimaculata, 


revealed by extensive field surveys over the last 
decade ( Ren et al., 2001 ; Qiu et al., 2004 , 2007b; 
Li et al., 2011). 

Among the three dominant Encarsia species, E. 


& Polaszek is a 


and heteronomous autoparasitoid. 


bimaculata Heraty solitary , 
arrhenotokous , 
Mated females lay eggs internally in B. tabaci 
nymphs, and develop as primary parasitoids. Males 
develop as hyperparasitoids, either in females of 
their own species or in other primary aphelinid 
parasitoids ( Antony et al., 2004). This parasitoid is 
widely distributed in China, Australia, India, 
Thailand, Philippines, Sudan and the USA ( Heraty 
and Polaszek, 2000). E. formosa is a cosmopolitan 
species ( Polaszek et al., 1992) and reproduces by 
microbe-induced thelytoky and has been successfully 
control of Trialeurodes 


used in biological 


vaporariorum on greenhouse crops in many countries 
(van Lenteren et al., 1997; Hoddle et al., 1998). 
E. sophia is also a cosmopolitan species found 
worldwide and a dominant parasitoid of B. tabaci 
( Kajita et al., 1992; McAuslane et al., 1993; 
Stansly et al., 1997; Monti et al., 2005) and T. 
vaporariorum (Gerling, 1983). 

It is assumed that, when choosing natural 
enemy species as biological control agents, closely 
related species usually share the same habits and 
host preferences, and phylogenetic relationships 
might be used to assess desirable candidates 
(Babcock et al., 2001). However, so far very little 
is known about these parasitoids in China, mainly 
because Encarsia populations in China have been 
rarely included in such studies, which might have 
delayed the accurate assessment and utilization of 
Encarsia parasitoids in whitefly control. 

In the present study, the 28S rRNA D2 and D3 
expansion regions of Chinese E. bimaculata, E. 
formosa and E. sophia populations were sequenced 
and analyzed with sequence data of other populations 
around the world to elucidate the relationships among 
them. The 


results are considered valuable in 
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evaluating the potential of local sources of these 
Encarsia parasitoids according to their relationship to 
the populations already reported, when choosing 
natural enemies in biological control programs. 


2 MATERIALS AND METHODS 


2.1  Parasitoid specimens 

The study was largely based on the materials 
collected in 2006 -— 2009 from the North, East, 
South and Southwest of China as well as from Kuala 
Lumpur, Malaysia (2007), and Egypt (2008 - 
2009). All samples of the parasitoids were collected 
from the nymphs of whitefly in the fields, and each 
sample was noted with collecting date, locality and 
collector. The nymphs of B. tabaci parasitized by 
parasitoids were kept in emergence chambers and the 
emerged parasitoids were captured and transferred to 
95% ethanol for storage. The whitefly species were 
identified either based on pupal morphology or by 
comparing their mitochondrial cytochrome oxidase I 
(mtCOI) sequence with that reported for B. tabaci 
(Qiu et al., 2007b). The emerged parasitoids were 
identified by Dr. Huang Jian (Fujian Agriculture 
and Forestry University, China). In addition to Æ. 
bimaculata, E. sophia and E. formosa, several other 
species of Encarsia and Eretmocerus that emerged 
from the B. tabaci nymphs in our surveys were also 
included in the study. All voucher specimens of 
Encarsia and Eretmocerus are deposited in the Insect 
Museum of Department of Entomology, South China 
Agricultural University, Guangzhou, China. 
2.2 DNA extraction and sequencing 

The DNA extraction was done according to the 
description of De Barro and Driver (1997) with 
double 


individual 


minor modifications. 
distilled water to 

parasitoids were homogenized in 200 uL lysis buffer 
(1% SDS, 10 mmol/L Tris-HCl, pH 8.0, 25 
mmol/L EDTA, 25 mmol/L NaCl, proteinase K 200 


mg/mL) using a 0.5 mL micro centrifuge tube and 


After washing by 


remove alcohol, 


micro pestle. The homogenate was incubated at 56°C 
for 2 -3 h in a water bath and then at 95°C for 10 
min to inactivate the Proteinase K. After incubation, 
samples were centrifuged for 1 min and directly used 
for PCR amplification or stored at -20 for 
later use. 

Sequences of the 28S rRNA D2 and D3 
expansion regions were also obtained from 14 
populations of 4 Encarsia species, 2 Eretmocerus 
species as well as from the outgroup Amitus 
hesperidum ( Table 1), three individuals of each 
selected for sequences 


species were randomly 
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analysis. The PCR primers and amplification 
programs of 28S rRNA D2 and D3 expansion regions 
(2000). All 
PCRs were performed in a 25 uL reaction volume 
and included 2.5 mmol/L MgCl, , 200 mmol/L of 
each dNTP, 1 pmol/L of each primer, 1 unit DNA 


Taq polymerase (Invitrogen, Guangzhou, China). 


were described in De Barro et al. 


Samples were amplified using a thermocycler ( Bioer 
XP cycler, TC-XP-D). Afterwards, 5 uL of PCR 
products were electrophoresed respectively using 
1.0% agarose gels with 1 x TAE at 8 V/cm for 2 h, 
then the gels were dyed in 10 mg/mL ethidium 
bromide for 30 min. When bands with expected size 
were visible in the gels, the other 20 uL of PCR 


products were used for sequencing. 


2.3 Data analysis 

Twenty five D2 and ten D3 sequences of 
Encarsia parasitoids were downloaded from GenBank 
and used as reference sequences (Table 2), and 
analyzed together with the sequences generated by 
us. All DNA sequences of 28S rRNA D2 and D3 
regions of aphelinid parasitoids were first analyzed 
and aligned with DNAStar (Lasergene® v5.0) and 
ClustalX 1. 83 ( Thompson et al., 1997). The 
phylogenetic relationships of different parasitoid 
populations were performed using NJ ( Neighbor- 
joining) options available in MEGA 4. 0. Genetic 
distances based on the D2 sequences were calculated 


with the Kimura 2-parameter model of MEGA 4.0. 


Table 1 Parasitoid species of Bemisia tabaci used in the phylogenetic analysis 


Parasitoid species 


Encarsia bimaculata 

E. bimaculata 

E. bimaculata 

Encarsia formosa 

E. formosa 

Encarsia japonica 

E. japonica 

Encarsia sophia 

E. sophia 

E. sophia 

Eretmocerus mundus 

E. mundus 

Eretmocerus sp. nr. furuhashii 
Eretmocerus sp. nr. furuhashii 
Amitus hesperidum 

E. bimaculata 

E. bimaculata 

E. bimaculata 

E. formosa 

E. formosa 

E. japonica 

E. japonica 

E. sophia 

E. sophia 

E. sophia 
Eretmocerus mundus 


Eretmocerus sp. nr. furuhashii 


Collecting locality 


Guangzhou, China 
Ismailia, Egypt 
Kuala Lumpur, Malaysia 
Kunming, China 
Hengshui, China 
Fuzhou, China 
Guangzhou, China 
Guangzhou, China 
Ismailia, Egypt 
Kunming, China 
Cairo, Egypt 
Ismailia, Egypt 
Guangzhou, China 
Sanya, China 
Guangzhou, China 
Guangzhou, China 
Ismailia, Egypt 
Kuala Lumpur, Malaysia 
Kunming, China 
Hengshui, China 
Fuzhou, China 
Guangzhou, China 
Guangzhou, China 
Ismailia, Egypt 
Kunming, China 
Ismailia, Egypt 


Guangzhou, China 


28S rDNA 
GenBank accession no. 
D2 D3 
+ JF819982 
+ JF819983 
+ JF819984 
+ JF819985 
+ JF819986 
+ JF819987 
+ JF819988 
+ JF819989 
+ JF819990 
+ JF819991 
+ 
+ JF819992 
+ 
+ 
+ JF819993 
+ JF819994 
+ JF819995 
+ JF819996 
+ JF819997 
+ JF819998 
+ JF819999 
+ JF820000 
+ JF820001 
+ JF820002 
+ JF820003 
+ JF820004 
+ JF820005 


+ : Target DNA sequences of the related parasitoid used in the phylogentic analysis. The same below. 
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Table 2 Parasitoid species used as references in the phylogenetic analysis 


Parasitoid species 


Encarsia bimaculata 
bimaculata 
bimaculata 
bimaculata 
bimaculata 
bimaculata 
bimaculata 
bimaculata 
bimaculata 


bimaculata 


bimaculata 
Encarsia formosa 
. formosa 

. formosa 

. formosa 

. formosa 

. formosa 


. formosa 
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. formosa 

Encarsia guadeloupae 
Encarsia sophia 

E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E. sophia 
E 


. sophia 


3 RESULTS 


Collecting locality 


Indial 
Taiwan, China 
Australial 
Australia2 
Sudan 
Guatemala 
India2 
Israel 

USA 
Australia 
Australia 
Egypt 
Greece 
UK1 

UK2 

USAI 





USA2 
Australia 
USA3 
Tonga 
Hawaii 
India 
Pakistan 
Ethiopia 
Florida, USA 
Spainl 
Malaysia 
Thailand 
Spain2 
Australial 
Australia2 
Tahiti 
Ethiopia 
Australia 


Polynesia 


3.1 Phylogeny of 28S rRNA D2 expansion 
segment 


The amplified region of the 28S rRNA D2 


expansion gene was about 590 bp in all of the 


28S rDNA 
D2 


D3 


GenBank accession no. 


AF254214 
AF254215 
AF254216 
AF254217 
AF254218 
AF254219 
AF254220 
AF254221 
AY599395 
AY615748 
AY615749 
AF223372 
AF223373 
AF223374 
AF223375 
AF223377 
AY599387 
AY615759 
AY615760 
AY615762 
AF254195 
AF254196 
AF254197 
AF254199 
AF254200 
AF254201 
AF254202 
AF254203 
AF254204 
AF254205 
AF254206 
AF254207 
AY599393 
AY615743 
AY615744 
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aphelinid parasitoids and sequences from the 2 — 3 


replicated individuals of the same species were 


identical to each other. 


Because of the short 


fragment in some of the reference sequences retrieved 


from GenBank, all D2 sequences were trimmed after 


alignment to a final length of 459 bp. 


The phylogenetic 


tree of these 


aphelinid 
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parasitoids based on 28S rRNA D2 expansion region 
is shown in Fig. 1. In the E. sophia clade, the 
populations from Hawaii ( GenBank accession no.: 
AF254195 ) and Tahiti ( GenBank accession no.: 
AF254207 ) were clustered into one subclade firstly , 
then they were clustered into one higher subclade 
with all other E. sophia populations including two 
from China. In the E. bimaculata clade, all the E. 
bimaculata populations except one from Australia 
(GenBank accession no.: AF254217) were clustered 
into one subclade. After that, the E. sophia clade 
and the E. bimaculata clade were clustered into one 
clade, which is a sister clade of E. japonica, the 
latter also included two E. japonica populations from 
Fuzhou and Guangzhou, China, respectively. 
Finally, the E. sophia, E. bimaculata and E. 
japonica formed a sister clade of the E. formosa 
clade. The phylogenetic tree also showed that E. 
sophia and E. bimaculata are sister species to each 
other, and they are more closely related to E. 
japonica than to E. formosa. 
3.2 Phylogeny of 28S rRNA D3 expansion 
segment 

The amplified region of the 28S rRNA D3 
expansion gene was about 312 bp in all of the 
aphelinid parasitoids, and the same length of 
sequences were aligned together with some of the 
GenBank. The 
phylogenetic tree of these aphelinid parasitoids based 
on 28S rRNA D3 expansion region is shown in Fig. 


reference sequences from 


2. In the NJ tree, three E. bimaculata populations 
from the USA ( GenBank accession no.: 
AY599395 ), Australia ( GenBank accession no.: 
AY615749 ) and Kuala Lumpur, Malaysia were 
firstly clustered together, followed with another 
Australian population ( GenBank accession no.: 
AY615748 ). 
were clustered with the E. bimaculata population 
China, then with the £. 
population from Egypt. 
Meanwhile, in the E. sophia clade, two populations 


Afterwards, these four populations 


from Guangzhou, 
bimaculata Ismailia, 
from Africa, the Ethiopia population ( GenBank 
accession no.: AY599393 ) and the population from 
Ismailia, Egypt were first clustered together, while 
the other four populations from Australia (GenBank 
accession no.: AY615743 ), Polynesia ( GenBank 
accession no.: AY615744) and from Kunming and 
Guangzhou, China were clustered closely. One 
population of Encarsia guadeloupae from Tonga 
(GenBank accession no.; AY615762) was clustered 
as one sister subclade of five populations including 
two from China ( Hengshui and Kunming). The 
phylogenetic patterns among the interspecies of 


Encarsia parasitoids based on 28S rRNA D3 gene 
were consistent with that based on 28S rRNA D2 
gene in Fig. 1. 
3.3  Intra-species genetic distances of the three 
Encarsia species 

The intraspecific genetic distances of the three 
Encarsia species were calculated using the 28S rRNA 
Of all the 459 positions in the 


alignment, 11 positions varied, which were located 


D2 sequences. 


between positions 16 and 437. The intra-species 
genetic distance of E. sophia is shown in Table 3. 
There are two base pair differences between the 
Guangzhou and Kunming populations of E. sophia 
(positions 173 and 235) , both populations have the 
same genetic distance with the other E. sophia 
populations. For example, they have the closest 
genetic distance with the populations of Ismailia 
Egypt, Australia 1 and 2, Spain 1 and 2 as well as 
the population from Ethiopia (0. 0022), then they 
have a relatively further genetic distance with those 
populations from India, Pakistan, Malaysia and 
Florida USA (0.0044 ), followed with larger 
distance with Hawaii and ‘Tahiti populations 
(0.0066). The largest genetic distance was found 
between the E. sophia CGuangzhou/ Kunming 
population and the Thailand population (0. 0088 ). 
In E. bimaculata, 8 positions varied in all the 
459 positions of the 28S rRNA D2 region, and all of 
them were located between positions 70 and 369. 
The sequences of the Guangzhou, Ismailia Egypt, 
Australia 1, Sudan and Guatemala populations are 
identical to each other, thus they have the same 


other E. 


being closer to the populations of 


genetic distance to the bimaculata 
populations, 
Taiwan, India 1 and 2, Israel, Kuala Lumpur and 
Australia 1 ( 0.0022 ) than to the Australia 2 
population ( 0.0066 ). Within E. formosa, the 
populations of Hengshui and Kunming from China, 
UK, USA (GenBank accession no. AF223377) and 


Greece are identical to each other, and have the 
same genetic distance to the population from Egypt 


(0.0022). 


4 DISCUSSION 


Several extensive surveys have reported 41 
species of Encarsia parasitoids of whitefly pests in 
China (Ren et al., 2001; Qiu et al., 2004, 2007b; 
Li et al., 2011). Although numerous studies have 
been reported on the morphological and/or molecular 
identification of Encarsia parasitoids ( Polaszek et al., 
1992, 1999; Schauff and Evans, 1996; Huang and 
Polaszek, 1998; De Barro et al., 2000; Babcock et 
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En. sophia Pakistan (AF254197) 
En. sophia Malaysia (AF254202) 
En. sophia Guangzhou China 

En. sophia Florida USA (AF254200) 
En. sophia Australia (AF254206) 
En. sophia Thailand (AF254203) 
En. sophia Hawaii USA (AF254195) 
En. sophia Tahiti (AF254207) 

En. sophia Kunming China 

En. sophia India (AF254196) 

En. sophia India (AF254196) 

En. sophia Spain (AF254204) 

En. sophia Spain (AF254201) 

En. sophia Australia (AF254205) 

En. sophia Ethiopia (AF254199) 
bimaculata Australia (AF254217) 
bimaculata Taiwan China (AF254215) 
bimaculata Australia (AF254216) 


bimaculata Ismailia Egypt 


En. bimaculata Sudan( AF254218) 


En. 


bimaculata Guangzhou China 


En. bimaculata Guatemala (AF254219) 


En. 
En. 
En. 


En. 


bimaculata India (AF254214 ) 
bimaculata Israel( AF254221 ) 
bimaculata India (AF254220 ) 


bimaculata Kuala Lumpur Malaysia 


100_| En. japonica Fuzhou China 


En. japonica Guangzhou China 


En. formosa Hengshui China 


En. formosa Kunming China 


100 


En. formosa UK (AF223375 ) 


En. formosa USA (AF223377 ) 


En. formosa Greece (AF223373 ) 


En. formosa Egypt (AF223372 ) 


En. formosa UK (AF223374 ) 


100 


Er. sp. nr. furuhashii Guangzhou China 


99 Er. sp. nr. furuhashii Sanya China 


100 | Ær. mundus Australia (AF273667) 


91 || Er. mundus Ismailia Egypt 


Er. mundus Cairo Egypt 


Amitus hesperidum Guangzhou China 


The Neighbor-joining tree constructed with Kimura 2-parameter model showing the phylogenetic relationship of 
different aphelinid parasitoids based on their 28S rRNA D2 expansion region sequences 
Values above the branches indicate the percent of 1 000 bootstraps. The same for Fig. 2. 
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En. bimaculata USA (AY599395) 


66 


En. bimaculata Australia (AY615749) 


62 


En. bimaculata Kuala Lumpur Malaysia 


69 


En. bimaculata Australia (AY 615748) 


61 


En. bimaculata Guangzhou China 


62 En. bimaculata Ismailia Egypt 


64 


85 


64 
70 


En. sophia Ethiopia(A Y599393 ) 
En. sophia Ismailia Egypt 


En. sophia Kunming China 


En. sophia Guangzhou China 


En. sophia Australia (A Y 615743 ) 


En. sophia Polynesia (AY615744) 


99 | En. japonica Fuzhou China 


En. japonica Guangzhou China 


95 


En. guadeloupae Tonga(AY615762 ) 


En. formosa USA (AY599387) 
99 En. formosa Australia(A Y615759) 
En. formosa Kunming China 


En. formosa USA (AY615760 ) 


En. formosa Hengshui China 


99 Er. sp. nr. furuhashii Guangzhou China 


Er. mundus Ismailia Egypt 


l 
0.005 


Fig. 2 The Neighbor-joining tree constructed with Kimura 2-parameter model showing the 


phylogenetic relationship of different aphelinid parasitoids based on their 28S rRNA D3 expansion region sequences 


al., 2001; Schmidt et al., 2001; Pedata and Polaszek , 
2003; Schmidt et al., 2006; Qiu et al., 2006, 
2007b), few of them examined the phylogenetic 
relationship between the Chinese Encarsia populations 
with populations from other countries. This has 
prevented us from getting a relatively accurate 
evaluation of the biology of either the Chinese or 
introduced Encarsia populations without conducting 
extensive investigations. 

In the current study, the intra-species phylogeny 
of Encarsia parasitoids based on 28S rRNA D2 and D3 
expansion regions are consistent with each other, 
phylogenetic of E. 


bimaculata, E. sophia, E. formosa as well as E. 


revealing the same pattern 


japonica. The neighbor-joining tree based on D2 
sequences revealed a range of genetic divergence of 
0 -17.88% , while the D3 phylogenetic tree revealed 
a range of genetic divergence of 0 -4.64% , indicating 
that the 28S rRNA D3 expansion region is more 
conservative than the D2 region. Our result is 
consistent with the report of Schmidt et al. (2006 ) 
that D3 is the most conversed in the 28S rRNA D2, D3 


and ITS-1 genes. Therefore, we subsequently used in 


this study the 28S rRNA D2 gene segment to 
investigate the intraspecific genetic distance of the 
three dominant Encarsia species in China. 

Our results on the intra-species genetic distance 
analysis of the three dominant Encarsia species 
indicated that there is genetic divergence between 
different populations within a single species, and the 
genetic 


distance is not consistent with their 


geographical ranges. For example, (Guangzhou 
population of E. sophia is closer to those populations 
from Australia, Spain, Egypt and Ethiopia, but further 
to the populations from India, Pakistan and Thailand. 
Also, Guangzhou population of E. bimaculata is closer 
to those populations from Egypt and Sudan, but further 
to populations from Taiwan. Furthermore, even in 
Australia, two populations of E. bimaculata (GenBank 
accession no.: AF254216, AF254217) have different 
genetic distance from other E. bimaculata populations. 
The worldwide these 
Encarsia populations through international trade and 
transportation of ornamentals with parasitized whitefly 


nymphs on the undersurface of plant leaves may be the 


distribution of commercial 


reasons for the inconsistency in the genetic distance 
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and geographical range of aphelinid parasitoids. 

To summarize, our current result of genetic 
divergence within each species is more informative for 
evaluating the biology of these biocontrol agents. Our 
results also indicated that, compared to the tree of 28S 
rRNA D3 which has more resolution for the intra- 
specific relationship study of Encarsia parasitoids, the 
phylogenetic tree of 28S rRNA D2 region is less 
informative and not of much value to evaluate the 
biological characteristics of closely related aphelinid 
parasitoids, however, the D2 tree can be more useful 
in species identification of aphelinids since it provided 
highly robust support for each species. Furthermore, 
the reference data of 28S rRNA D3 region of Encarsia 
in GenBank are limited in numbers so far, more 
populations should be sampled for the systematic and 
comprehensive investigations that focus on the intra- 
and interspecific relationships of Encarsia wasps in 
future studies. 
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